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Abstract 

In this paper we provide a prescription for solving the initial value 
constraints of GR for algebraically general spacetimes in the polar 
representation. The prescription uses the Ashtekar magnetic field and 
two eigenvalues of the CDJ matrix as inputs, and provides as an output 
the SO(3,C) frame solving the Gauss' law constraint, in the form of a 
triplet of SO(3,C) angles. The two eigenvalues are the physical degrees 
of freedom, which fix the third eigenvalue through the diffeomorphism 
and Hamiltonian constraints. We provide two fixed point iteration 
procedures for determining the angles, as well as an analysis of the 
algebraically special cases. 
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1 Introduction 



The Gauss' law constraint of general relativity G a is a triple of differential 
equations signifying the invariance under internal SO(3,C) rotations. In 
Paper VII we have applied this constraint to the rectangular representation 
of the CDJ matrix which serves as the momentum space variable of the in- 
stanton representation, thus reducing it by three degrees of freedom. In this 
paper we chose the diagonal (anisotropy) elements as the physical degrees 
of freedom, and used the Gauss' law constraint to establish the shear (off- 
diagonal symmetric) elements as images of the anisotropy elements under 
the action of the Gauss' law propagator Jj, seen as a map. It was found in 
the rectangular represenation that the Hamiltonian constraint amounts to 
a nonlinear relation amongst these degrees of freedom requiring an iterative 
procedure in order to disentangle. But in the polar representation, it turns 
out that the eigenvalues of the CDJ matrix more suitably characterize the 
physical degrees of freedom of the instanton representation when one con- 
siders the Hamiltonian constraint. So in the present paper we will in similar 
fashion use the Gauss' law constraint to establish a map from the eigen- 
values to the corresponding nonphysical degrees of freedom. These latter 
degrees of freedom are the parameters which define the SO(3, C) frame. In 
both representations the propagator map requires some input based on the 
choice of coordinates used for the inversion. This latter input is directly en- 
coded in the choice of connection Af, which is more fundamentally derived 
from a triple of congruences of linearly independent integral curves 7 which 
fill 3-space S. The tangent vector fields to these curves determine B l a , which 
can in turn be used to find Af (see e.g. Paper VI). 

The outline of this paper is as follows. In section 2 we write the initial 
value constraints in polar form starting from the Ashtekar variables. This 
identifies Xf = (Ai,A2,A2), the eigenvalues of ^( ae \, as the physical degrees 
of freedom. There are three unknowns, namely a triplet of SO(3, C) angles 
9, which define a 50(3, C) frame where ^> ae is annihilated by G a . In section 
3 we provide a prescription for determining the angles 9 to first and second 
order using the Gauss' law constraint combined with equivalence of the 
rectangular and polar representations of ^ ae . The degeneracy of eigenvalues 
determines which angles are fixed by this prescription, which angles are 
undetermined, and also the algebraic classification of the CDJ matrix. We 
have associated this algebraic classification to the Petrov classification of the 
corresponding spacetime. For algebraically general spacetimes there is no 
obstruction to computing the angles to first order. In section 4 we provide 
a fixed point iteration procedure to determine the angles to all orders, when 
a solution exists. The solutions are labelled by two eigenvalues (Ai, A2) and 
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7, the congruence of integral curves in £ whose tangent vectors v a = B l a di 
are constructed from the Ashtekar magnetic field B l a . Section 5 provides an 
alternate iteration procedure by phrasing the Gauss' law constraint directly 
as a set of differential equations on the angles 9. 
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2 Initial value constraints problem in polar form 



The initial value constraints of general relativity take on a polynomial form 
in the Ashtekar variables [3], [5]. The diffeomorphism constraint is given by 

Hi = e ijk P a B k a = 0, (1) 

which signifies the invariance under spatial diffeomorphisms in S. The 
Hamiltonian constraint signifies invariance under deformations normal to 
£ and is given by 

H = ±e ljk e abc a a o{d k c + e ijk e abc a^ b B k = 0, (2) 

where A is the cosmological constant. The Gauss' law constraint, which 
signifies invariance under left-handed SU(2)_ rotations on internal indices, 
is given by 

G a = D t al = d t al + f abc A b a\ = 0, (3) 

where f a j, c are the structure constants for the left-handed SU(2)- algebra. 1 
The CDJ Ansatz [7] 

K = *aeBl (4) 

was introduced by Riccardo Capovilla, John Dell and Ted Jacobson as a 
means for solving (1) and (2). The CDJ matrix f ae G SO(3, C) <g> SO{3, C) 
can be parametrized by its antisymmetric part a ae and its symmetric part 
A ae , which can in turn be parametrized by a polar decomposition, 2 

iv = a ae + \ ae = e aed ?p d + o af {e)\ f oJ e (e). (5) 

Here, A = (Ai, A2, A3) are the eigenvalues of \ a b = *(ae) an d tpd is a SU(2)- 
valued 3-vector, while O ae € SO(3,C) implements a complex orthogonal 
transformation parametrized by three complex angles 6 = (#i,#2>#3)- in 
exponential form this is given by O = e d ' T , where T are 50(3) generators 
satisfying the Lie algebra 

1 For notational purposes, we will frequently refer to the Gauss' law constraint as G a 
and to the Hamiltonian and diffeomorphism constraints as H M = (H,Hi). 

2 This requires the existence of three linearly independent eigenvectors so that the 
symmetric part ^( a b) is diagonalizable. 
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[Ta,Tb] = ie A BcTc- 
Under (4), the diffeomorphism constraint (1) yields 



(6) 



Hi = e i]k BiB k e ^ ae = e ljk BiB k e (e aed ^ d + A ae ) = V x, (7) 

which implies that ip d = 0, or that CDJ matrix is symmetric = ^^ a b)- 
For nondegenerate B a (7) reduces to 

e ljk B{B k e e aed ij d = (detB^B' 1 )^ = 0, (8) 

which implies that ip d = 0. 

Substitution of (4) reduces (2) to 

H = (dct5) {^Var^ + Adettf ) = Vx (9) 

where we have defined Var^f = (tr\I/) 2 — tr^ 2 . Substitution of (5) into (9) 
upon dividing by det-B ^ 0, which requires nondegeneracy of B l a , yields 

H = War(X ae ) + Adet(A ae ) + (AA ae - 5 ae )M e = 0. (10) 

The first two terms of (10) can be rewritten explicitly in terms of the eigen- 
values of \ ae due to the cyclic property of the trace, yielding 

H = (AiA 2 + A 2 A 3 + A 3 Ai) + AA1A2A3 

+ (AO af O ef \ f -5 ae )ip a ip e = 0. (11) 

Therefore, on the space of solutions to the diffeomorphism constraint for 
nondegenerate B l a , the terms in (11) quadratic in ip d vanish and all ref- 
erence to the 50(3, C) angles 9 become eliminated from the Hamiltonian 
constraint. Equation (10) then simplifies to the following algebraic relation 
amongst the eigenvalues A 

a + ^ + ^ + ^ = ^ a ' = -( aa iA , + a 1 + J - (12) 

Hence for Aj 7^ (12) fixes one eigenvalue A3 completely in terms of the 
remaining two eigenvalues Ai and A 2 , which is independent of B a . 3 

3 We will ascribe to Ai and A2 the interpretation of the physical degrees of freedom of 
the momentum part of the phase space of GR. 
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While the 50(3, C) angles 9 have been eliminated from the Hamiltonian 
constraint on solutions to Hi = 0, we will see that they appear explicitly in 
the Gauss' law constraint G a . Therefore G a will be regarded as a condition 
using eigenvalues A/ satisfying (12) and the connection Af as inputs to fix 
9, which defines an equivalence class of SO(3, C) frames. The eigenvalues 
A f must then be rotated from the intrinsic frame where 9 = into the frame 
where ^ ae becomes annihilated by the Gauss' law constraint. This frame, 
fixed by the correctly chosen 9, should correspond to the solution of the 
initial value constraints. 

Upon substitution of (4) into (3) we have 

Dial = VaeDiBt + B* A*oe = ^A*ae, (13) 

where we have used the Bianchi identity DiB l e = 0. This leaves remaining 
the action of the covariant derivative on the CD J matrix, which takes its 
values in the tensor representation of SU(2)- 

A*ae = di^ ae + fabcA^ce + * ac^/e&c- (14) 

Defining vector fields v e = B\di and the 'helicity density matrix' Cb e = 
A\B l e , the Gauss' law constraint then reduces to 

W e {* ae } = B*A*oe = V e {1V} + {fabfSge + febgS a f)C be ^f g = 0. (15) 

Unlike which are algebraic, G a are a set of differential equations which 
can be written in the form 

w e {\ f (e- 6 - T ) fa (e- 6 - T ) fe } = 0. (16) 

Prom (7) and (12) we have obtained a two parameter family of solutions 
to = (H, Hi) per point labelled by Ai and A2 in the polar representation 
of * ae . Recall that O = e 9 ' T , where T are the generators and 9 = (9 1 ,9 2 , 9 3 ) 
are a set of complex rotation parameters. Hence given 

/ Ai \ 

Kb = O ae {9) A 2 Oj b (9) 
V A 3 J ef 

where Ai, A2 and A3 are related by (12), we must implement the Gauss' 
law constraint. Equation (16) is a set of three coupled nonlinear first order 
differential equations for the three angles 9, one equation for each index 
a. The solution to these differential equations is an outstanding problem, 
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which constitutes the last remaining obstacle to the construction of a general 
solution to the initial value constraints of general relativity [8], [9]. In this 
paper we will provide a prescription for obtaining the 50(3, C) angles 6. 
The CDJ matrix can be parametrized in the rectangular form 



*ae = e f ae <p f + E f ae V f , (17) 

f f 

where e ae and E J ae form a basis of diagonal and off-diagonal symmetric 3 by 
3 matrices. We have shown in Paper VII using the Cartesian representation 
that the solution of G a establishes a map (p — > ^ [<p\ from the anisotropy to 
the shear elements, reducing ty ae by three degrees of freedom. In this basis 
the Gauss' law constraint is given by 



which can be written explicitly in matrix form 
v i — C[23] —C32 C23 

C3I v 2 — C[31] —C13 

— C21 C\2 V3 — Cr 12 l 





0. 






II 






K<P3 J 



(18) 



(V2 + C13 — 2C31 C22 — C33 V3 + 2C21 — C\2 \ / *3 \ 

Vl + 2C 32 - C 23 v 3 + C21 - 2Ci 2 C 33 - C u *i = 

Cn - C 22 v 2 + 2C 13 - C31 Vl + C7 32 - 2C 23 / \ *2 / 

where C[ j,] = (7^ — Cfe a . In paper VII it is shown how the above operator- 
valued matrix equation can be formally inverted to yield 



*/ = Jffg = -(^w e ) 4 e>^ 9 ). (19) 

Jj defines the propagator from the anisotropy to the shear subspace, using 
the former as a source term and the connection Af as an input. Hence 
iff can be regarded as the physical degrees of freedom analogously to the 
eigenvalues Aj. For integrable configurations Af the propagator Jj is well- 
defined. For other configurations, Jj is path-dependent. Irrespective of the 
issue of integrability, part of the results of this paper will be to provide the 
map <p — > A along with its physical interpretation, which should correspond 
to a solution to the Gauss' law constraint. 4 

4 For certain configurations this solution is well-defined and for other configurations it 
is not well-defined. 
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2.1 Physical interpretation of the £0(3, C) frame 

Define 7 as the congruences of the integral curves of three vector fields vi , 
V2 and V3, which are linearly independent at each point i£E. The vector 
fields v a define three internal directions at each point, and are independent 
of coordinates as is 7^ , and we will regard 7 as the fundamental structure 
of 3-space S. Choose coordinates x l = (x 1 , x 2 , x 3 ), and define the three by 
three matrix B % a , given by 



K = VaW}. (20) 

The columns of B l a , seen as a three by three matrix, define three directions 
in 3-space S. Associate the spatial index i to vector components in a spatial 
manifold S, and associate the index a to vector components in an internal 
6*0(3,(7) manifold. Let us place indices i and a on the same footing and 
split B ia into its symmetric and its antisymmetric parts 





(* 


c 






B a i — 


C 


y 








V b 


a 








The symmetric part £( a i) = \{B a i + Bi a ) defines a triple of column 3- vectors 
s a with components i, labelled by the internal index a. Seen as a triple of row 
vectors in internal space with components a and labelled by the spatial 
index i, B % a defines the same triple 













: 


| = S'l. s 2 = 1 




| = S 2 > s 3 = 1 


: 






V a j 







si = ^ c j = Si, s 2 = ^ y J = s 2 , s 3 = ^ a j = s 3 . 

Hence, the axes defined by B^ with respect to the 3-manifolds E and 
50(3,(7) are the same. There are six degrees of freedom in 5( oi ), three 
D.O.F. fix the relative scaling of the axes, and the other three D.O.F. fix 
their orientation. 

The antisymmetric part B^ = \{B a i — Bi a ) also determines a triple 
of vectors, seen either as spatial column vectors a a labelled by the internal 
index a or as internal row vectors a!, labelled by the spatial index i 





ai = \ -w , a 2 = u , a 3 



which contains three degrees of freedom, and 



a' 1 = 




















l.-l-l 








V u ) 




T 
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Whereas for there is no distinction between s and s', for B^ the 

vectors a and a' differ by a reflection. So an arbitrary three by three matrix 
B l a defines two coordinate systems, s = (s!,s 2 ,S3) and a = ±(a!,a 2 ,a 3 ). 
The effect of applying a SO(3, C) transformation to B\ 

B'\ = (e e - T ) ae Bl, (21) 

is to rotate the internal index e. If B l a is taken to be a magnetic field then 
this rotation, which is a SO(3, C) gauge transformation, rotates s and a by 
an amount encoding the existence of a SO(3, C) frame labelled by 6. 5 



5 The identification of B' a as a magnetic field is equivalent to the existence of a connec- 
tion A1 for which B l a satisfies the Bianchi identity DiB l a = 0. 
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3 General solution for the 50(3, C) angles 



We have shown on the diffeomorphism constraint shell that the S0(3, C) 
angles 9 become eliminated from the Hamiltonian constraint H which then 
depends only on the eigenvalues A of *( ae ). But it is clear that 9 explicitly 
appear in the Gauss' constraint G a . Our interpretation is to view G a not as a 
condition on A, but rather a condition that fixes 9 once the physical degrees 
of freedom have been specified. We will regard the eigenvalues A as the 
physical degrees of freedom, in a similar way that <p are the physical D.O.F. 
in relation to in the Cartesian representation of *S> ae in (18). Therefore, to 
have a solution to the initial value constraints in the polar representation, 
one may to establish a map <p — > A from the anisotropy elements of ^ ae 
to the eigenvalues. As we will see, this process when convergent fixes the 
angles 9 which in turn identifies an equivalence class of 50(3, C) frames 
corresponding to each Af. 

The angles 9 may be interpreted as the parameters of a 50(3, C) trans- 
formation of the diagonal matrix of eigenvalues Diag(\\, A2, A3) into a new 
frame. This is given by 

*ae = {e e - T ) af \ f {e^ T ) fe . (22) 
Expanding out the commutators in (22), we have 



*ae = SaeK + [6- T, 5 ae X e ] + ±[6- T, [9 ■ T, 5 ae X e }} + .... (23) 

where (T A ) ae = e Aae are the generators of 50(3). The zeroth order term of 
(23) is given by X ae = 8 ae X e . The linear term is given by 



[9- T, 8 ae X e ] = 9 A ({T A ) ae \e - X a (T A ) ae ) = 9 A (T A ) ae {\ e - A ). (24) 
The second order term of (23) is given by 



[9-T,S ae Xe 



±9*9 A 



(TB)af(T A ) fe (X e - A/) - (T B ) a f(Xf - X a ){T A ) fe 



^9 B 9 A (T B ) af {T A ) fe (\ e - 2X f + A ). (25) 



The n th order term is given by 
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^e A ^ A ^...e A "(T Al ) aai (T A2 ) aia2 

n 

• • • (TA n - 1 )a n - 2 a„-i{TA„)a n - 1 a n <)a„e — l^Cfc A afc , (26) 



fc=0 



where = n\/k\(n — k)\. To prove this we will use mathematical induction. 
First make the following definition for the product of generators 

(TA 1 A 2 ...A n )aa n = (TAi )aoi (Ta 2 )a 1 a 2 • • • (Ta (27) 

In order to compute the n+l th term we will need to evaluate the commutator 
[6 A T A ,T Al A2...A n \- Note from (26) that the n th term can be written as 

TA 1 A 2 ...A n - 1 A n (CqX^ - C?A a2 + C*2 A ffl3 - . . . (-l) k C^\ ak+1 + . . . (-l)"C , ™_ 1 A an ^ 

= T A2 A 3 ...A n A n+1 (Cq \ a2 - C"A a3 + C% X a4 — ... {-l) k C^\ ak+2 + ■■■ (— l) n C'"_ 1 A an+1 ^ , (28) 

where we have relabelled the dummy indices to start from Ai vice A±. To 
obtain the commutator, subtract the left-multiplication of the the right hand 
side of (28) by (T Al ) aai from the right-multiplication of the left hand side 
by (T An+1 ) an a n+1 , obtaining 

9 A T A , T Al A 2 ...A n (Co\ ai - C"A a2 + C% A a3 - . . . (-l) fc C£ A afe+1 + . . . (-l) n C'™_ 1 A an ^ 

= (^41)001(2^2)0102(^43)0203 • • • (^A„)a„_ia„(?i„ + i)(i„a„ + i 

(Cq A ai — C"A a2 + C% X a3 — C3 A a4 + . . . (— l) fc C^A afc+2 + . . . (— l) n C"_ 1 A a+ i^ 

~~ (Tai)ooi (^4 2 )aia2 (^4 3 )o 2 o 3 • • • (^4„)a n _ia„ ( T A n+1 ) 
(Co^a 2 — C™A a3 + C2 A a4 — • • • (-l) k C%\ ak+2 + . . . (-l) n C"_ 1 A a+ i^ 
ra n+1 

= (T Al A 2 ...A n+1 ) £(-1)*^ + = ^M 2 ...A„ +1 ^(-l) fc C fc " +1 A afc (29) 

fc=0 fc=0 

where we have defined = for k < and have used the identity C^_ x + 
C™ = C™ +1 . Contraction with # Al # A2 . . . An 9 A ^ and dividing by (n + 1)! 
provides the desired result. Hence we have a prescription for evaluating the 
expansion to any order desired. The full expansion of the CDJ matrix is 
given by 

^ae = S ae \ e + A (T A ) ae (\ e - X a ) 
00 1 N N N 

+ e n ^ n (^ja^o™ ^(-i)"cao„, (30) 

N=2 ' k=l m=l n=0 
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where the remainder term contains terms of second order and higher. 

To establish the map from the shear elements to the angles we 
must equate the polar representation of ^f ae to its Cartesian representation. 
Hence we must solve the relation 

*ae = (e - T )abXb(e- 9 - T ) be = e^Vf + E^ f . (31) 

The Gauss' law constraint in the Cartesian representation given by (18), 
with solution (19) which implies the following form for ^ ae 

* ae = (el + E[ e J g f )^ = T° e <p g , (32) 

where Jj = —(M~ v )f a H ag is the Gauss' law propagator acting on the di- 
agonal CDJ matrix elements <pf. We will show in this paper that subject 
to (32), equation (31) saturates the degrees of freedom in 9 A , which corre- 
sponds to the choice of a SO(3, C) frame where the constraints are satisfied. 
The polar representation of the CDJ matrix \I> ae can then be seen as a rota- 
tion from the diagonal frame of eigenvalues into this SO(3,C) frame where 
it becomes annihilated by the Gauss' law constraint. We would like to ex- 
press 9 A as a functional of Ai and A2, as well as the configuration chosen for 
the Ashtekar connection Af. Equation (31) may be written 

Ai 9 3 (X 2 -X 1 ) 9 2 (X 1 -X 3 ) \ 
3 (A 2 -Ai) A 2 9 1 (X 3 -X 2 ) +...= 
9 2 (X 1 -X 3 ) 9\X 3 -X 2 ) A 3 / 

j§<Pa (02 Jlfa 
J$tp a Jf(f a 993 / 

where the dots denote all terms of second order and higher. The zeroth 
order term is independent of the angles 6 A = 9^, which implies ^ = X a 
for a = 1,2,3. The result is that to zeroth order, one may equate the 
diagonal CDJ matrix elements to the eigenvalues. The following relation 
ensues to first order 

(A 3 - x 2 )ol 1) = J{x r , (x 1 -x 3 )ef 1) = jix r ,(x 2 -x 1 )ef 1) = j(x f . (33) 

When the eigenvalues are unequal Ai ^ X 2 ^ A3, then one may divide (33) 
by the eigenvalue differences. Then the shear CDJ matrix elements, which 
are symmetric, provide the first nontrivial approximation to the angles A 
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9 »> = ^(^) J ^- <34) 

where I a b c = \e a j, c \ is an even version of the epsilon symbol. Before proceed- 
ing with the recursion relation for 9 let us first treat the situations where 
the eigenvalues are equal, such that (34) cannot be written from (33). 

3.1 Algebraic classification of the CD J matrix 

One way to organize the solutions to the initial value constraints is accord- 
ing the algebraic classification of the CD J matrix \l/ ae . Knowledge of the 
eigenvalues of the Weyl curvature tensor for a spacetime determines the 
range of possibilities for its algebraic classification, up to knowledge of its 
eigenvectors. Define the CDJ matrix by 

Ke=~Sae + ^ae, (35) 

where ipae is symmetric and traceless in its indices. Now define a triple of 
objects T\\ B , which form an isomorphism between internal internal SO(3, C) 
indices a = (1,2,3) and unprimed SL(2,C) spinorial index pairs AB = 
(00,01, 11). Then the traceless part of can be written in the form 

</W = Va B Ve D i>ABCD, (36) 

where ipABCD = 'P(abcd)-, which is totally symmetric in its indices. If we 
take ipABCD to be the spinorial decomposition of Weyl, the antiself-dual 
part of the Weyl curvature tensor, then ip ae contains the right number of 
degrees of freedom to represent Weyl in SO(3, C) language. When diago- 
nalizable, the left and right hand sides of (35) can be diagonalized by the 
same 50(3, C) transformation 

(e / )oe(^) =^ae + A / / (e^) ae (37) 

with no summation over /. In (37) Aj are the eigenvalues of ip ae and A/ are 
the eigenvalues of Vtv. The association of ip ae to Weyl then implies that 
the algebraic classification of spacetime is also encoded in ^> ae in SO(3,C) 
language. The Hamiltonian constraint can be solved point by point in S to 
furnish the eigenvalues Ai, A2 and A3 of ^ a ei of which according to (12) two 
independent ones Ai and A2 are freely specifiable. From (37) it is clear that 
the eigenvalues of ipABCD and of are directly related. The multiplicity 
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of the eigenvalues of ipABCD is also freely specifiable which in turn deter- 
mines the Petrov type within a given class, up to the number of linearly 
independent eigenvectors. So by solving the Hamiltonian constraint, one 
narrows down the possibilities to (I, D, O), (II, N) or III . 

However, more information is needed in order to uniquely fix the space- 
time. In addition to the eigenvalues one must also know the number of lin- 
early independent eigenvectors, then one can construct a solution. For any 
triple of eigenvalues Xf there exist an infinite number of 50(3, C) frames 
(22) to transform into. However, only certain frames correspond to a GR 
solution, which is the purpose of the Gauss' law constraint to determine. 
One the angles 6 are found, then Weyl be explicitly constructed using (36) 
and (35). It is then a simple matter to deduce the algebraic type of the 
spacetime by finding the eigenvectors of ^> ae . 

3.2 Algebraically special case with one distinct eigenvalue of 
the CDJ matrix 

In this case all three eigenvalues equal Ai = A2 = A3. This corresponds 
to spacetimes of type O (three distinct eigenvectors), type N (two dis- 
tinct eigenvectors) and type III (one distinct eigenvector). According to 
the Hamiltonian constraint, 

£ + A_ -A,~f. (38) 
Since all three eigenvalues are equal, the eigenvalue differences are given by 

Ai - A 2 = A 2 - A 3 = A 3 - Ai = 0. (39) 

According to the general iterative procedure for Ai = A 2 = A3, the angles to 
first order must satisfy 

WAc-A 6 )0f 1} = J f A X f . (40) 

But the left hand side of (40) vanishes on account of the equality of the 
eigenvalues, which means that 9^ cannot be determined using the iteration 
procedure. The right hand side of (40) is given by 

j s f <p g = (E b fa w b )- 1 (e c ag w c )<p g . (41) 

Note that when one sets ip g = X g , as required for the zeroth order of iteration, 
that {e c ag w c )(p g vanishes. This can clearly be seen in matrix form from the 
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anisotropy contribution to the Gauss' law constraint upon setting <p g = X g = 
Ai for g = 1,2,3 



v l — C[23] —C32 C*23 

C31 V2 — C[ 31 ] — C13 

— C21 C12 V3 — C[!2] 



Ai > 




^ V1A1 


Ai 


1- 


v 2 Ai 


Ai J 




\ v 3Al 









Note, due to the definition C< ae \ = C ae — C ea that the contributions due to 
the magnetic helicity density matrix C\, e cancel out, leaving v a A = B l a di\ = 
Va. Multiplication by yields Ai = for i = 1, 2, 3 since i?* is non- 

degenerate, or that the single eigenvalue Ai to linearized order is spatially 
homogeneous. Hence based on G a alone, Ai = A2 = A3 = *S>(t), correspond- 
ing to arbitrary 9. Then ^ ae is the isotropic CDJ matrix ^f ae = S ae ^>, where 

= *&(t) is an arbitrary function of time, independent of position in S. 
Consistency with the Hamiltonian constraint then requires that = — -| 
as in (38), which is a numerical constant. 

It follows that 6^ is undetermined by the Gauss' law constraint in this 

case. Since the angles 9 are undetermined at linearized level, the 50(3, C) 
rotation matrix O ae (9) G 50(3, C) is also undetermined. This indetermi- 
nacy should be preserved to all orders of the iteration procedure which we 
will describe, implying that all 50(3, C) frames provide a solution, and 
are equivalent. This yields three sets of possibilities forming the equiva- 
lence classes of spacetimes for this case. A nondegenerate O ae contains of 
three linearly independent eigenvectors, which corresponds to Petrov type 
O. When O ae contains either two or one linearly independent eigenvectors, 
one should expect spacetimes of types iV or /// respectively. 

3.3 Algebraically special case with two distinct eigenvalues 
of the CDJ matrix 

In this case Ai = A2 / A3. According to the Hamiltonian constraint, we 
have 



2 1. 

— + — + A = 0^A 3 
Ai A3 

The eigenvalue differences are given by 



■f Al V 

VAiA + 2/ 



(42) 



A3 — Ai — A3 — A2 



-Ai 



AiA + 3 



,AiA + 2 

Therefore, for the angles 9 at linearized level we have 



(43) 
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0i _ 1 f AiA + 2\ 



9 3 
^(1) 



? 2 



A + 3, 



Arbitrary; 
Ai 



1 / AiA + 2 n 
Ai VAiA + 3/ 



(A1A + 2) 



(4 + ^)Ai 



f3/ 
J 2 



(44) 



One angle is undetermined to first order. While there are two distinct 
eigenvalues in this case, all three eigenvalues are fixed by Ai through (42). 

To first order, # 3 is arbitrary and (9 1 ,9 2 ) are determined by Ai. Upon 
application of the iteration to the next and to all orders, the angles 9 should 
have the following functional dependence 



= 01%^,^]. (45) 

Evidently, the angles depend upon two degrees of freedom not counting the 
connection, which implies that \l/ ae = \l/ ae [0^ ; Ai , Af] also inherits those 
degrees of freedom. 

Since B l a is freely specifiable, then one should in general have A3 ^ 
Ai, which restricts the number of linearly independent eigenvectors to 2 
or 3. All freedom has been exhausted with the exception of # 3 , which is 
freely specifiable. A choice of 9 3 corresponding to three linearly independent 
eigenvectors should then take on the interpretation of type D spacetimes, 
while two linearly independent eigenvectors should yield type II. 



3.4 Algebraically general spacetimes with three distinct eigen- 
values 

Algebraically general spacetimes, spacetimes of type I, have three distinct 
eigenvalues Ai 7^ A2 7^ A3 for the CDJ matrix ty ae and correspondingly three 
distinct eigenvectors. The Hamiltonian constraint for this case is given by 

r + r + f + A = ' ( 46 ) 

Ai A2 A3 

which furnishes A3 function of Ai and A2 

A3 = "(aa 1 a 2 + A a 1 + a 2 )- (47) 

The eigenvalue differences are given by 
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/AAi + 2 + ^\ /AA 2 + 2+^ , 

A3 - Al ^(^Ft) ;A3 - A2 ^(^Ft'' (48) 

and the 50(3, C) angles to first order are given by 




J3A1 + Jf A2 — Jf 




(49) 



All three angles can be found to first order, hence for Type I spacetimes 
it should in principle be possible to continue the iteration to all orders, for 
each configuration 



3.5 Iteration to second order for the algebraically general 
case 

For the algebraically general result of Ai / A2 7^ A3, the next step 

is to iterate the solution for 6 to the next order. To linearized order we have 

We will now evaluate the second order solution, incorporating the remainder 
term to this order. Recall from (25) that this is given by 

= ^9 B 9 A (T B ) af {T A ) fe (\ e - 2X f + A ). (51) 

We need to compute the elements of the matrix corresponding (51) to sec- 
ond order. An efficient way is to draw a lattice, with all indices arranged 
along the horizontal and a duplicate set along the vertical axis. Then one 

6 The configuration B l a are derived objects from v a upon the choice of a coordinate 
system x l , via the relations B % a = v a {i'}. The vector fields v a are tangent to the set 
of (linearly independent) triples of congruences of integral curves 7, which fill 3-space 
E. Since we regard 7 as fundamental structure, then the 50(3, C) angles are functionate 
9 = 0[Ai, A2; 7] of the eigenvalues and the integral curves 7. 
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tabulates the trivial and nontrivial contributions due to epsilon symbols for 
each possible pair, for a given starting symmetry. 

Using this scheme, one deduces from (51) that for the diagonal terms 
a = e, only the B = A terms contribute while the B ^ A terms vanish. 
Likewise, for the off-diagonal terms a / e one finds that it is now the A ^ B 
terms which contribute with the A = B terms vanishing. Hence, an easy 
way to visualize the interaction amongst the indices is that (A, B) mimic 
the symmetries of (a, e). 7 Writing out the terms for explictness, one has for 
the diagonal terms that 



a = e = 3 — ► 9 1 9 1 (X 2 - A 3 ) + e 2 9 2 {\ 1 - A3); 
a = e = 2 — > 3 3 (Ai - A 2 ) + 9 1 9 1 (X 3 - A 2 ); 

a = e = 1 — ► e 2 e\x 3 - Ai) + e 3 e 3 (\ 2 - (52) 



and for the off-diagonal terms that 



a = 1; e = 2 — > 0^ 2 (Ai + A 2 - 2A 3 ); 
a = 2; e = 3 — > ^ 3 (A 2 + A 3 - 2Ai); 
a = 3; e = 1 — > 9 3 9 1 (X 3 + Ai - 2A 2 ). (53) 



We will need to solve the equation 

MA e + e efg 9 f 9 9 (X g - A e )) + I fae 9 f (X e - X a ) 
+ ^Ifae0 a O e (X a + A e - 2X f ) = el eVf + E&frd- (54) 

Hence, one has for the zeroth order term that 

v4 0) = A/- (55) 
Then, one substitutes (55) to find the shear elements, which gives 9^ 



^> = 7 ^(vh;)« A /- (56) 

We must now compute the new anisotropy elements to this order 

^ = MAe + eefgd^/^Xg - A e )). (57) 



7 This observation is useful, particularly when computing the higher order terms by 
hand. 
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Let us show a sample calculation for illustrative purposes 



•P? = Ai + (l^-) 2 (A 3 - A t ) jf A„ - (^-L^) 2 (A, - A 2 ) jjA„ 

- (*• + (*T^r)* " (aT^)*)*- (58) 

Note that there has been a cancellation of one factor of the eigenvalue dif- 
ferences from the denominator. In shorthand notation, the general term is 
given by 



5fe ~ tdgf 



1 



^A a — A y 

Moving on to the error terms, we have 



J e f )X, 



(59) 



,(i) 



m - IabcQ(i)0(i)(^b + Ac - 2A.4). 



(60) 



Putting in the expression from (56), this yields 



p(i) _ ( A& + A c — 2A a 

* a -2!^ c V(A c -A a )(A a -A fe ) 



(61) 



Hence, to obtain the second order solution we must solve the equation 



1 



2!^ fec V(A c -A a )(A a -A 6 ) 



Afe + A c — 2A a 



f,9 



The angles are given by 



to first order, and to second order by 



(62) 



(63) 



labc 



Ab + A c — 2A a 



(A c - A a )(A a — Ab) 



(64) 



For integrable configurations A" the process can be repeated iteratively to 
arbitrarily high order, which should lead to a solution to the initial value 
constraints of GR to the chosen order. 
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4 Fixed point iteration for the 50(3, C) angles 



We will now derive a recursion relation to find the angles to any desired 
order, starting from eigenvalues Aj and the Ashtekar connection A?. 8 The 
projections into anisotropy and shear subspaces are given, respectively, by 
(ff = e a ^ae and = EJ e ^ a e- At the same time, tyf = ^ f[<p] constitutes 
a map between these subspaces under the Gauss' law constraint, which we 
require to be preserved at each order of the iteration. Let the full solution 
to all orders be given by 



* ae = A ae [A, B a ] = (^ A >> T ) a6 A 6 (e-^> T ) be , (65) 

which is a functional of A/ and A\. Then (65) must be annihilated by the 
Gauss' law constraint in the sense of (18) and (19), such that 



tt, = Ef^ ae \\ A?} = e« e J|^ ae [A, A?} = Jfa g . (66) 

The effect of this is to establish a map iff — > ^ f[<f\- 

Hence, we set up a recursion as follows. The CDJ matrix can be written, 
at the level prior to implementation of the Gauss' law constraint, as 



)ae 



* ae = (e d - T Xe~ e - T ) c 
= e Aae e A (X e - X a ) + E[ e Ef'(e e - T Xe- e - T U 
-e A a e e A (\e ~ Xa) + e{ e ef (e^Ae'^W (67) 

with summation over a'e' but not over ae. At this point all we have done is 
to decompose the CDJ matrix into its orthogonal subspaces in a Cartesian 
basis and to add and subtract eAaeQ A {X e — A a ) = to the equation. But 
note that this isolates the term linear in 9, which will provide a linearized 
term for the expansion. Next, write the Gauss' law constraint ^/ = J fag 
in the following form 

Ef{e 9 - T Xe- 9 - T ) a>e> = ef'jy- T Xe- 9 - T ) a , e , (68) 
Project (68) into the shear subspace by contraction with e{ 



E{ e Ef\e^Xe-^) a , e , = ef E{ e Jf {e^ Xe^ T ) a , e , (69) 

8 The main principle is the observation that the CDJ matrix elements corresponding 
to the zeroth and the first order term in the polar representation fall into orthogonal 
subspaces, which enables one to define a fixed point iteration procedure with a convenient 
physical interpretation. 
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The left hand side of (69) is the same as the first three terms on the right 
hand side of (67). Hence, we can write down the relation 

-e A aeO A (X e ~ K) = E{ e ef Jj (e^ \e- e - T ) a , e , . (70) 

We have now written the angles 6 A explicitly in terms of themselves, and 
therein lies the basis for a recursion relation. Transfer the second and third 
terms of (70) to the right hand side 

£AaeO A {K — A a ) = eAaed A (K ~ A a ) 

+ (^^-^)(W T W. (71) 

Assign the index n + 1 to the left hand side of (71) and the index of n to 
the right hand side and we obtain, upon dividing by £Aae(A e — A a ) ^ 0, 9 

C-i) = fa +lA ™(x^x) E ^ E 1 e ' ~< e 'j 9 f)(^y T \e-^y T ) a , e , (72) 

The physical interpretation is that the righmost term of (72) in brackets 
corresponds to the CDJ matrix ^ ae = at the n th order of iteration 

y{n) _ ^ [n yT Xe -e {ny T )a ^ (73) 

(n) 

While ^ae automatically solves the Hamiltonian constraint H for all n due 
to choice of eigenvalues according to (12), it does not solve the Gauss' law 
constraint G a for n < oo. In this limit we would have 

hm^^f' - ef'ty^l = 0, (74) 

which implies that the operator E^ e — ej e J| in its action on provides a 
measure of how far away one is from the exact solution. This error is then 
applied to the angles 6*( n ) to produce the next level of refinement 6*( n+ i)- 
Then is used to recompute the corrected CDJ matrix *ie +1 ^ using 

the same eigenvalues A which solve the Hamiltonian constraint H by con- 
struction. The process is then repeated in cycles ad-infinitum. 

As a double-check on this result, if one chooses #( ) = 0, then (72) yields 

9 That this is possible arises from the fact that the zeroth and the first order terms in 
the expansion of \& ae [A;6] fall into orthogonal subspaces. 
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6 (i) = lA ™{x^) E ™ E f e '(<'e' X ^ ~ lAae( J ^)ELe*' e 'j 9 f (e h a , e ,\ h ).(75) 

Note that for d^L = 0, the second term on the right hand side of (75) vanishes 
on account of the orthogonality of the shear and anisotropy subspaces 



£&£7? V (<WV) = 0, (76) 

and (72) reduces to (34). Hence, the process has produced the correct result 
to first order and can then be extended to all orders. Hence, the angles in 
the algebraically general case are given by 9 A = \\m. n ^ > . oo A ^y where = 
and the recursion relation (72) holds. 

Hence, we regard the Gauss' law constraint as being a constraint on 9 
but not a constraint on Xf. The eigenvalues Xf form a direct input in the 
determination of the angles 9 a = 9 a [X; B], which defines a special 5*0(3, C) 
frame where all the constraints are satisfied. The iteration procedure then 
proceeds in cycles using two operations. First one rotates the eigenvalues 
Xf) into the 50(3, C) frame for the given order, followed by computation of 
the resulting error and adjustment to the corrected 50(3, C) frame for that 
order. Then one repeats the process ad-infinitum. 

4.1 Necessary and sufficient condition for the existence of a 
fixed point 

Arguments for convergence of a fixed point for the recursion (72) can be 
motivated as follows. The norm of satisfies the following bound 

|*/,(r)| < |( e T - e[Al ' A2;7l )//'||(e T - 0[Al ' A2;7l ) flfl '4 fl ,||A^| < e 2 ^'|A| (77) 
Vf e S, where 



|A| = VlAil 2 + |A 2 | 2 + |A 3 | 2 ; |0| = vW + |0 2 l 2 + |0 3 l 2 (78) 

When the 9 are purely real, then O ae = (e ) ae is a unitary transformation 
fixing \^f g \ to |A|, the norm of the eigenvalues. In this case the iteration 
procedure must converge. When 9 is complex, then the existence of a so- 
lution is directly related to the sign of the real part of 9 = Re[9] + ilm[9]. 
For Re[9] < the iteration map is a contraction, which should yield a fixed 
point. For Re[9] > 0, additional steps are required to establish the existence 
of a fixed point, if one exists. 
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We now present an argument for the convergence of the SO (3, C) angles 
under (72) as follows. Recall that the CDJ matrix * ae is given by 



* ae = vPaelAi.Aa;^ = (e^- A2 ^- T ) a/ (e^- A2 ^- T ) e/ A / . (79) 



On the left hand side of (79) is the CDJ matrix solving the constraints, 
which is explicitly a functional of two eigenvalues Ai and A2 and the integral 
curves 7 of v a , and on the right hand side is the polar parametrization of 
the same matrix. On the left hand side is a Cartesian representation of * ae 
with respect to three Cartesian directions in an internal 577(2) _ <g> 5C/(2)_ 
space labelled by indices a, e. The ability to write down this Cartesian 
representation implies the finiteness of the right hand side of (79), since the 
right hand side has been obtained by the multiplication of three matrices. 

If one could compute * ae by an alternate method, then one could es- 
tablish finiteness of the elements on the left hand side of (79). This would 
imply finiteness of the right hand side, which in turn would imply that 6 
are well-defined and hence the fixed point iteration used to determine them. 
Since Diag(\i, A2, A3) is the diagonalized form of * a e, then the full matrix 
can be constructed from its eigenvectors if there are three linearly indepen- 
dent eigenvectors. The CDJ matrix, which is symmetric, can be written in 
the form 



The characteristic equation for * ae is given by the cubic polynomial 



which determines the eigenvalues X p ~ (Ai, A2, A3). The solution for to (80) 
is given (see appendix D for the derivation of the solution to the cubic) by 




det(* ae - (5 ae A) = A 3 - (ipi + <p 2 + y?3)A 2 - (<pi<P2 + ¥>2¥>3 
- *i - *1 - *i)A " (2*1*2*3 - <Pi<P2P3) = 0, 



(80) 




r* / 4 ( Vary (tgg \ p 



4(det* + 



(tv^){Var^) _ 2(trfr) 3 ^' 
3 27 / 



^A ^Var^ + (tTg)3 ^ 



(81) 



where p = 1, 2, 3 



label the three roots, which serve as the eigenvalues A, and 



T 1 / 3 (x) = cos [-cosx]. 



(82) 
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However, rather than use the CDJ matrix to find its eigenvalues we will 
do the reverse. Starting from the eigenvalues A, we will construct the CDJ 
matrix ^> ae from scratch. 

The eigenvalues of the CDJ matrix are given by the equation X ae u e = 
\5 ae u e , where u e are the eigenvectors. This can be written as 



(fiu + + = \ p u; 
^&3n + (f2V + f iro = X p v; 

^2U + i&iv + 1P3W = XpW, (83) 
which leads to the solution, labelled by p, of 

1 \ 

^2^3+^i(A p -yi) 

*l*3+*2(Ap-^2) 

(Ap-yi)(Ap-y 2 )-1 f 3 

X ' P \ *l*3+*2(A p -¥. 2 ) / p 

for some normalization factor u p . Recalling the polar parametrization of the 
CDJ matrix 




*ae = O ab 0){e{ c X f )O T ce {e) (84) 

and the fact that a symmetric matrix can be diagonalized by the orthogonal 
matrix consisting of a set of linearly independent eigenvectors, we then have 
the interpretation of the rotation angles as in 

Vi v 2 v 3 J . 
W\ W2 w 3 J 

The ability to write down a solution for the eigenvectors then implies con- 
verge of the fixed point iteration procedure. 10 



4.2 Construction of a solution 

We have presented a prescription for determining the angles 6 to all orders. 
Equation (72) refines the estimate of the angles 9 A using a combination of 
6*0(3,(7) transformations and Gauss' law inversions. The fixed point, if it 
exists, is then given by 9 A = lim^oo^^ = ^[Ai, A2; 7]. 

To convert the quantities computed into an explicit solution, one must 
first reconstruct the CDJ matrix using 6 

10 The application is the following. One freely chooses three distinct eigenvalues Ai, A2 
and A3 solving the Hamiltonian constraint for each point r £ E and from this determines 
the 50(3, C) rotation matrix simply by construction of the corresponding eigenvectors. 
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*oe = *ae[Ai,A 2 ;7] = (e^^- T ) o/ A / ( e -^ A ^- T ) /e (85) 
where the eigenvalues A/ are related via 

A3 = "(aa 1 a 2 + A \ 1 + a 2 )- (86) 

Equation (85) and (86) consitute a general solution to the constraints parametrized 
by the two physical degrees of freedom (Ai, A 2 ). n From ^ ae [A, B a ] one may 
determine the radiation properties of the spacetime by construction of the 
Weyl scalars via 



*ae = " + 1>ae — > */ = Sf tt" 1 • (87) 

where S^ 6 is a projection operator onto the five components of a massive spin 
two tensor representation. ip ae is the left-handed part of the Weyl curvature 
tensor, expressed in SU(2)_ language. Vae is symmetric and traceless, and 
can be expressed in the language of self-dual two forms 



IpABCD = ^aeVAB^CD- ( 88 ) 

Additionally, one may reconstruct the metric tensor by first defining the 
internal metric rj ae from the CDJ matrix solving the constraints, as in 



Vae = WAi, A 2 ; 7] = (det^)(*- 1 ^- 1 ) ae . (89) 

Equation (89) contains information regarding the SO(3,C) frame, which 
solves the Gauss' law constraint G a . Additionally, it contains information 
regarding the choice of coordinate system used to solve G a . The 3-metric is 
obtained by projection of (89) into 3-space S via 



MAi,A 2 ;7] = VaeiB-^UB-'TjidetB). (90) 

The matrix B l a , which encodes the choice of coordinate system, may be seen 
as a gauge group element such that hij is simply r] ae in a different gauge. 
Defining the one form 

11 From Ai and A 2 one can determine the algebraic classification of the spacetime as 
alluded to earlier. 
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J = dx l + N l dt, 



(91) 



where = (N, N l ) are the lapse function and shift vector which are freely 
specifiable, one may reconstruct the spacetime metric 

ds 2 = g^daTdx" = -N 2 dt 2 + h^J ® J. (92) 
In terms of the internal variables, this is given by 

ds 2 = -N 2 dt 2 + r/ ae [Ai, A 2 ; j]oo a <g> oj b , (93) 

where co a = (det.B) 1 / 2 ^- 1 )^. 

Hence the spacetime metric as well acquires the label of its algebraic 
classification as encoded in Ai and A 2 . The CD J matrix can now be expressed 
in spacetime via 

%l P „ = ^%(~S ae + ^(e-^^le-^^), (94) 

where = ^f ae F^ u (by the CD J Ansatz), are the components of a set 
of self-dual two forms and F£ v = F® U [A] is the curvature two form for the 
self-dual Ashtekar connection Af. 

Note that the Ashtekar curvature B % a occurs explicity in the solution 
(94), along with the eigenvalues Xf. The two momentum physical degrees 
of freedom of GR are encoded within Ai and X 2 . There must be two degrees 
of freedom's worth of configuration variables conjugate to (Ai,A2) which 
must be accounted for. However, the only configuration variables in (94) 
reside within the Ashtekar curvature Bf. Since B l a in general contains nine 
components, it follows that the information contained within seven of these 
components must be regarded as unphysical. 12 



12 This is reasonable, since B l a is determined through (20) from a choice of coordinate 
system. 
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5 Gauss' law constraint in the polar representa- 
tion: Revisited 



We have provided a prescription for determining a special 50(3, C) frame 
in two stages, namely by solving the Gauss' law constraint in the Cartesian 
representation, and then transforming the solution for ^> ae into the polar 
representation. It is also possible to bypass the Cartesian representation, 
writing the Gauss' law constraint as a set of differential equations directly on 
the angles 0. This approach provides some insight into the intrinsic nature 
of Gauss' law and gauge transformations in the instanton representation 
variables. We would like to solve the equation 

G a = v e {* ae } + O^* /ff = 0, (95) 

where v e = B\di are a set of vector fields constructed from the magnetic 
field B l e . We have defined 

C l 9 = (daffbge + $gefabf)Cbe, (96) 

where C ae = A^B l e is the local 'magnetic helicity density matrix' of the con- 
nection Af = Af [a], 13 which in (96) appears in a kind of tensor representa- 
tion with f a b c the structure constants of 50(3, C). Since ^> ae is symmetric, 
it can be written in a polar decomposition when there are three linearly 
independent eigenvectors [10] 14 

*ae = (e e - T ) af \ f (e- e - T ) fe , (97) 

where 6 = (0 1 ,^ 2 ,^ 3 ) are three complex angles and T = (T J 4) ae = tAae 
are the generators of 50(3) in the adjoint representation. The complex 
orthogonal matrices have the property that 

(e 9 - T f ab = {e e - T ) ba = (e-^U; det(e^) = 1. (98) 

First let us focus on the parts of (95) involving vector fields, namely v e {A ae }. 
Acting on (97) with v e and using the Liebniz rule, we obtain 

13 We use the term 'magnetic helicity matrix' as a gravitational analogue to a quantity 
appearing in magnetohydrodynamics h = J E d 3 xA ■ B, the magnetic helicity of a Mawxell 
potential A with magnetic field B = V x A. While the latter object is global scalar 
measuring the linking of magnetic field lines in E, the gravitational counterpart Cae = 
A1B\ is a 50(3, C) x SO(3, C)- valued matrix defined locally at each point. 

14 The interpretation is a Lorentz transformation where Im[8\ and Re[8] are the rotation 
and the boost parts of the transformation, respectively. 
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v e {* ae } = v e {{e e - T )af\f{e- e - T )f e } = M(e eT )af}\f(e- e - T )fe 

+(e e - T )a f M^Ke" e - T )fe + (e e - T ) af X f M(e- e - T )fe}. (99) 

Continuing with the expansion we have 



MO A Ke°- T )ar(T A ) ff X f (e-^) fe + (e'^^WA/} 

-v e {9 A }(e e - T ) af X f (T A ) ff ,(e- 9 - T ) f , e 
= (e e - T ) af (e- e - T ) fe M^f} + MO A Ke e - T )af>(e e - T ) ef (T A ) rf X f 

-v e {0 A }(e e - T U(e e - T ) ef (T A ) ff X f . (100) 

Relabelling indices / -H- /'on the last term of (100) we have 



(e e - T ) a/ (e- e - T ) /e v e {A / } + (e e - T ) af (e°- T ) f/eVe {0 A }(TA)MX f - X f ) 



= {e e - T ) af v' f {X f } + (T A ) flf (X f -X f y f ,{e A } , (101) 



where the primed vector field is defined by 15 



= {e^ 1 ) de v e . (102) 

We now multiply (101) by (e~ ' T )h a , in anticipation of some further manip- 
ulations. This eliminates the prefactor and changes the index / — >■ h on the 
remaining terms. The result is 



(e- 6 - T ) ha Ve{Xae} = v'{A ft } + (T A ) f , h {X h - X f/ )V f ,{9 A }. (103) 
The second term in the last line of (103) decomposes into the following terms 

(T 1 ) 23 (A 3 - A 2 )v' 2 {0 1 } + (T 1 ) 32 (A 2 - A3)v 3 {^}; 
(T 2 ) 3 i(Ai - A 3 )v 3 {# 2 } + (T 2 ) 13 (A 3 - Ai)vi{0 2 }; 
(T 3 ) 12 (A 2 - Ai)vi{6» 3 } + (T 3 ) 2 i(Ai - A 2 )v' 2 {# 3 }, (104) 

which can be written as a matrix operator, with rows labelled by h and 
columns labelled by A, acting on 6 arranged as a column vector 



15 The vector field v e = B\di transforms as a vector under SO(3, C) gauge transforma- 
tions assuming that B r a transforms as a vector. Hence the primed objects correspond to 
an SO(3, C) transformation of the corresponding unprimed counterparts, parametrized 
by the angles 9. Therefore it should be understood from the context that the primed 
quantities are implicitly 8 dependent. 
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/ (A3-A1K (Ai-A 2 )V 2 \ / e 1 \ 

-Q hA e A = - (a 2 - a 3 )v' 3 o (Ax - a 2 )v; e 2 . 

\ (A 2 - A 3 )v' 2 (A3 - X 1 )V 1 J \0 3 J 

We have used the property that {Tj{) ae = — (T J 4) ea = e^ae- There is also a 
contribution from the first term on the second line of (101), given in matrix 
form by 

/ v ; \ / a 1 

v'JAU = v' 2 A 2 
V v' 3 / \ A 3 

We are now ready to move on to the magnetic helicty contribution to (95). 
5.1 Magnetic helicity density contribution 

Next, we expand the contribution C[ g Xf g in the polar representation. This 
is given by 

Cp*f g = (fabfSge + feb 9 S af )C be (e d - T )f 9 (e d - T ) gd Xd 
= (fabf(e 6 - T ) fd (e e - T U + f ebg (e e - T )ad(e e - T ) gd )C be \ d . (105) 

In a similar fashion to v e {A ae } we multiply (105) by (e~ 9 ' T )ha, yielding 

{e- 6T ) ha Cl*\ fg = (He e -- T ) ha (e- e - T ) df f afb (e 9 - T ) ed 

+fe bg (e~ 9 - T ) ha (e 9 ' T )a d (e e - T ) gd ) C be X d - (106) 
Using (e~ 9 ' T )f ia (e e ' T )a d = $hd this changes the index d into h, which yields 

(e- e - T ) ha Ci 9 X fg = (-(e- eT ) ha (e- e - T ) df f afb C be (e e - T ) ed X d 

+feb g C be (e 9 ' T ) g hXh = -(fh dc C' cd X d + C'^Xh). (107) 
In the last line of (107) we have defined 

C cd = (e~ dT )cbC be (e 9 ' T ) ed ; = (e~ e ' T )h g f gbe C be , (108) 

which can be seen as transformations of the unprimed quantities into the 
same 50(3, C) frame as v'. We have used 

{e- e - T ) ha {e- e - T ) df (e e - T ) bc f afc = f hdb , (109) 
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which transformed the last term of (107) into 



(e 9 ' T )ha(e 9 ' T )dffafbCbe(e 9 ' T ) e dXd 

= fhdc(e~ 9 ' T )cbCbe(e 9 ' T )edXd = fhdcC' cd \d- (HO) 
Equation (107) splits into the following terms 



fi23C 32 \2 + /132C23A3 + C^jAi; 

/23lCl2^3 + /213C 31 Ai + C[ 2 ]A2; 

/312C21^2 + /32lCl2^2 + C[3]A3, (HI) 

where we have defined 



C[i] = C[23]; C[2] = C[3i]; C[3] = C[i2]- (112) 

Hence (111) in matrix form is given by (incorporating) the minus sign from 
(107)) 





( - c m 


"^32 


r" 

°23 


(e- e - T ) ha Cf9\ fg = T hg \ g = | 


C31 


r" 

^[2] 


°13 




V -C21 


Cl2 


°[3] 




Combining all terms, we have that 

{e- 9T ) ha G a = -Q hA e A + D hg X g = (113) 
where we have defined 

D hg = S hg V g + F hg . (114) 

Equation (113) can be written in the following matrix form 

(A 3 -Ai)v 3 (Ai-A 2 )v 2 

(A 2 - A 3 )v 3 (Ai - A 2 )v; 

(A 2 -A 3 )v 2 (A3-A1K 



/-if r^l 

L/ 32 u 23 

C3I V 2 ~~ C[2] _ Cl3 

~C21 ^12 v 3 ~~ C[3] 







II 


'SI 




V 3 / 
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In the case where the eigenvalues are unequal Ai 7^ A2 7^ A3 the matrix Q A 
for certain configurations A® might be invertible. One could then solve for 
the angles using 16 The formal solution is the recursion relation 

A = (Q-y h D hg \ g , (115) 

Equation (97) formally establishes a map from the eigenvalues A, for a given 
configuration Af, to the angles 9. The solution requires the inversion of 
QhA, which is a matrix whose elements consist of noncommuting differential 
operators. The point is that (95) is not a constraint on A, but rather it is a 
constraint on 9. One may define a recursion relation 

9^ +1 = jf[A;9 n ]X g , (116) 
where 9q = and we have defined Gauss' law propagator 

jf[A;9 n ] = {Q- l [A-X9]) Ah D hg [A;9]. (117) 

The physical interpretation is that one chooses an SO (3, C) frame to start 
the iteration procedure in (115), and then computes the vector fields v a and 
magnetic helicity density C ae in this frame. Then (115), seen as a linear 
equation in 9 relative to the chosen frame, provides a correction. Then v a 
and C ae are recomputed for the corrected frame and re-inserted on the right 
hand side. The process is repeated ad-infinitum. 

It appears that a more in-depth analysis would be required to establish 
the convergence of this iteration procedure, however one can in principle 
determine whether the angles are bounded in the following sense. Applying 
the Minkowski inequality to (115), we have 

|0| < |Q _1 ||DA|. (118) 

The idea is that (115) and (116) can at any stage in the iteration procedure 
be reagrded as a linear equation in 9 n+ i, with a source term depending on 9 n . 
The 9 n dependence can be absorbed into the definition of the configuration 
variables B l a as a gauge transformation, and the equation seen as a linear 
equation in an unknown angle 9 forming the left hand side. One always has 
an equation of this form since it is written anew upon each iteration. Then 
(118) provides a bound on this angle based upon the choice of configuration, 
encoded in v a and C ae , and upon the eigenvalues. For the former one needs 
to determine a norm for the propagator, which should prescribe conditions 
for the functions A/ it is acting on. 

16 This is actually a recursion relation since the angles appear on both sides. 
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6 Appendix A: Cartesian representation of the CDJ 
matrix 

First decompose X ae into a basis of diagonal (anisotropy) and off-diagonal 
(shear) components 

\ ae = Ef e S f + ef ae <p f . (119) 

In (119) we have defined the following basis of symmetric three by three 
matrices, using e ae = fifa$fe f° r the diagonal and 

E\ e = <5 a2 5 e 3 + 5 a3 5 e2 ; E 2 ae = <5 a3 5 e i + 5 a i5 e3 ; E\ e = 5 a i5 e2 + 5 a2 5 e i (120) 
for the off-diagonal matrices, which satisfy the orthogonality relations 

EfEi e = 25 9 f e{ e ef = 5f g - e{ e Ef = 0. (121) 

This is known as the Cartesian representation. Further defining H ae = 
eieWf and M ae = E[ e Wf, the Gauss law constraint reduces to 



M af S f + H af <p f = 0^S f = -{M- l ) fa H ag ^ g = J 9 f <p g , (122) 

where Jj<f g = —(M~ 1 )f a H ag . Assuming the invertibility of M a f, (122) 
establishes the shear elements Sf as a map of the anisotropy elements (/?/, 
which upon solution to the Gauss' law constraint should reduce \ ae by three 
D.O.F. 17 



7 Appendix B: Inversion of the operator Q 

The inversion of the operator Q proceeds along similar lines as in the Carte- 
sian representation. The matrix that we wish to invert is of the form 

/ z y 
Q = w x 
\ v u 

We must find an operator- valued matrix U ae = 5 a \a e + b a2 b e + 5 a sc e , such 
that 

17 The specifics for the inversion of M a f are treated in Paper VII. Our purpose here 
merely is to outline the map from the physical D.O.F. to the angles 9, which can formally 
be defined irrespective of the well-behavedness of 9. 
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a 2 w + a 3 v aiz + a 3 u aiy + a 2 x 
& 2 w + 6 3 v 61 z + 6 3 u biy + 6 2 x 
C2W + C3V C1Z + C3U ciy + c 2 x 



is diagonal. Setting the off-diagonal parts to zero we have 



a\z + 0,311 = 


- 


-> (13 = 


— aizu - 


aiy + a 2 x = 


— 


-> a 2 = 


-aiyx" 


biy + 6 2 x = 


- 


->6i = 


-6 2 xy~ 


6 2 w + 63V = 


— 


+ 63 = 


— 6 2 wv~ 


c 2 w + c 3 v = 


— 


-> c 2 = 


— C3VW _ 


ClZ + C3U = 


= - 


— > Cl = 


C3 uz 



(123) 



We can now solve for the diagonal elements 



a 2 w + 03V = — ai(yx~ 1 w + zu _1 v); 
biz + 6 3 u = -6 2 (xy _1 z + wv 4 u); 
ciy + c 2 x = -c 3 (uz _1 y + vw 4 x). 



(124) 



This leaves the elements (a±, 6 2 , C3) undetermined. But these elements cancel 
from both sides of the equation yielding 




yx _1 w + zu _1 v \ / 1 — yx _1 — zu" 1 

xy x z + wv x u — xy 1 1 — wv" 1 

uz _1 y + vw _1 x J \ —uz 1 — vw _1 1 

The second factor has the interesting property that the Mb a = 

which is the analogue for vector-field valued matrix of antisymmetry for 

c-numbered matrices. 



32 



References 

[1] T. Thieman 'On the solution of the initial value constraints for gen- 
eral relativity coupled to matter in terms of Ashtekar variables' Class. 
Quantum Grav. 10(1993)1907-1921 

[2] Riccardo Capovilla, John Dell and Ted Jacobson 'The initial value prob- 
lem in light of Ashtekar variables 

[3] Ahbay Ashtekar. 'New perspectives in canonical gravity', (Bibliopolis, 
Napoli, 1988). 

[4] Ahbay Ashtekar 'New Hamiltonian formulation of general relativity' 
Phys. Rev. D36(1987)1587 

[5] Ahbay Ashtekar 'New variables for classical and quantum gravity' Phys. 
Rev. Lett. Volume 57, (1986) 

[6] Newman E.T. and Rovelli, Phys. Rev. Lett. 69, 1300 (1992) 

[7] Richard Capovilla, Ted Jacobson, John Dell 'General Relativity with- 
out the Metric' Class. Quant. Grav. Vol 63, Number 21 (1989) 2325- 
2328 

[8] Ricardo Capovilla, John Dell, Ted Jacobsen and Lionel Mason 'Self- 
dual 2-forms and gravity' Class. Quantum Grav. 8 (1991)41-53 

[9] Riccardo Capovilla, John Dell, Ted Jacobson 'The initial value problem 
in light of Ashtekar 's variables' arXiv:gr-qc/9302020 

[10] Asher Peres 'Diagonalization of the Weyl tensor' Phys. Rev. D18, Num- 
ber 2 (1978) 



33 



